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We  develop  a  first-principles,  self-consistent  theory  of  the  far-infrared  (FIR)  electromagnetic 
response  for  electrons  confined  in  a  quantum  dot.  We  find  that  for  small  electron  number  the 
FIR  absorption  spectrum  corresponds  to  that  associated  with  parabolic  confinement,  i.e.,  absorp¬ 
tion  dominated  by  a  single  peak,  which  occurs  at  the  frequency  corresponding  to  the  interlevcl  , 
separation  of  the  parabolic  potential,  and  is  roughly  independent  of  For  large  electron  num¬ 
ber,  an  upward  shift  in  the  resonance  frequency  occurs  as  the  electron  density  probes  the  increas¬ 
ingly  nonparabolic  curvature  of  the  dot  potential.  Effects  of  an  applied  magnetic  field  arc  also  in¬ 
vestigated. 


With  recent  advances  in  nanofabrication  technology  it 
has  become  possible  to  confine  electrons  in  all  three  spa¬ 
tial  dimensions  in  semiconductor  structures  called  quan¬ 
tum  dots. 1  “4  Such  structures  are  analogous  to  atoms  but, 
in  place  of  an  atomic  potential,  electrons  see  the 
artificially  constructed  dot  potential.  Typically,  the  la¬ 
teral  electron  confinement  is  produced  either  from  period¬ 
ic  etching1-3  or  periodic  gating2-4  of  a  quasi-two- 
dimensional  (2D)  electron  gas.  In  the  first  case,  the  elec¬ 
trons  in  each  dot  are  confined  by  layers  of  positively 
charged  impurities  and,  possibly,  occupied  surface  states. 1 
In  the  second  case,  a  gate  voltage  can  effectively  be  viewed 
as  producing  positively  charged  disks  that  confine  elec¬ 
trons  to  the  dots. 

A  recent  experimental  study2  of  the  far-infrared  (FIR) 
electromagnetic  properties  of  quantum  dots  shows  that  for 
zero  magnetic  field,  the  FIR  absorption  spectrum  is 
governed  by  a  single  peak  whose  corresponding  frequency 
is  roughly  independent  of  the  number  of  electrons  ne  occu¬ 
pying  each  dot.  For  nonzero  magnetic  field,  applied  per¬ 
pendicular  to  the  dots,  a  splitting  of  the  resonance  occurs 
producing  two  peaks  that  are  also  relatively  insensitive  to 
changes  in  nr.  This  behavior  can  be  explained5-8  when 
the  potential  that  confines  electrons  in  a  dot  is  parabolic. 
However,  the  parabolicity  of  the  dot  potential  can  be  ex¬ 
pected  to  break  down  when  the  electron  density  extends  to 
the  edges  of  the  dot.  Therefore,  a  complete  theory  is 
needed  of  the  electromagnetic  response  of  a  quantum  dot 
that  employs  a  realistic,  self-consistently  calculated 
confining  potential.  In  this  Rapid  Communication  we 
present  such  a  complete  approach.  Specifically,  we  calcu¬ 
late  self-consistently  the  FIR  electromagnetic  response  of 
a  GaAs  quantum  dot,  containing  between  I  and  30  elec¬ 
trons,  without  and  with  an  applied  magnetic  field. 


In  current  dot  geometries  the  dots  are  well  separated  so 
that  no  overlap  between  the  electronic  wave  functions  in 
neighboring  dots  can  occur.  Furthermore,  it  has  been 
shown8  that  the  interdot  electromagnetic  interactions  pro¬ 
duce  a  negligble  shift  in  the  FIR  absorption  frequency  for 
most  samples  studied  so  far.  We  therefore  assume  that 
each  dot  is  isolated  from  all  others.  Also,  since  the  lateral 
extent  of  the  electron  confinement  in  the  dot  is  typically 
about  ten  times  larger  than  that  along  the  growth  direc¬ 
tion,  we  model  each  dot  as  a  two-dimensional  system.  We 
also  approximate  the  positive  charges9  to  be  uniformly 
distributed  into  a  jellium.  Similar  modeling  has  been  used 
in  the  case  of  small  metallic  clusters, 10  where  electrons 
were  taken  to  be  confined  by  a  spherical  jellium.  We  fur¬ 
ther  assume  that  the  dot  potential  possesses  a  circular 
symmetry.  This  latter  assumption  has  been  shown,  even 
for  square-shaped  dots,”  to  be  a  reasonable  approxima¬ 
tion.  Finally,  for  mesa-etched  dots,  electrons  that  may  be 
trapped  in  surface  states1  are  modeled  by  a  ring  of  nega¬ 
tive  charge  at  the  edge  of  the  jellium  disk. 

The  potential  energy  K(x)  of  an  electron  interacting 
with  a  charge  density  p(x)  in  the  dot  is  given  by 

K(x)  “*—  f  t ,  (1) 

kj  |x  x  | 

where  x  —  (r,r),  r  ”,(r,0),  and  k  is  the  dielectric  constant 
of  the  semiconductor.  For  a  jellium  disk  of  radius  R  and 
areal  density  n+  centered  at  the  origin  of  the  r  "0  plane 
the  charge  density  is  p+  xmen+S(:)6(R  —  r ).  For  a  ring  of 
N  -  negative  charges  in  this  plane  at  r-/?  the  correspond¬ 
ing  density  is  p-  —  —  (eN  -/2xr)8(r— R)§(:).  The  po¬ 
tential  energy  of  an  electron  interacting  with  these  respec¬ 
tive  charge  densities  in  the  r  “0  plane  ts 
K(r,0)  —  V+{r)  +  V -(r),  where 
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In  the  equations,  K  and  E  are  the  complete  elliptic  in¬ 
tegrals  of  the  first  and  second  kind,  respectively. 

As  a  consequence  of  the  circular  symmetry,  the  wave 
functions  for  electrons  in  the  z  -0  plane  confined  by  the 
dot  potential  defined  by  Eqs.  (2)  and  (3),  and  including 
a  magnetic  field,  B*-Z?z,  are  separable:  Vn.i(r) 
” with  y/„,/  satisfying  the  Schrodinger  equa¬ 
tion: 


h '  1  d  d  ,  ,  . 

2^7  7 T/Tr  +K"<r) 


Petr (r)wmV(r)-\ - 7  +  t m,.(o}r~  +  r  Ifico,  ,  (4) 

2  m,r~ 

/“0,  ±  1,  ± 2. ,  0),. — —  . 

mec 

Here,  me  is  the  electron  effective  mass  and  /  is  the  orbital 
quantum  number.  The  Zeeman  splitting  of  the  electron 
levels  is  ignored.  The  potential,  K(r)“F+(r)  +V-(r) 
+  V,(r),  where  Ve(r)  is  obtained  by  solving  Eq.  (I) 


employing  the  electronic  charge  density  in  the  dot. 
p,  (r,z)  “  —  en(r)8(z),  with 

«5) 

«./ 

where  the  sum  is  over  all  occupied  levels  and  I,  or 

2  gives  the  occupancy  of  level  («,/);  we  take  the  tempera¬ 
ture,  7’“0  K.  The  ground  state  of  the  system  is  obtained 
by  self-consistently  solving  Eqs.  (I),  (4),  and  (5). 

We  now  calculate  the  charge  density  induced  in  the  sys¬ 
tem  by  external  electromagnetic  radiation  of  the  form 
K^'fr./)  —  KcU(r)exp(/ft>/),  where  a>  is  the  frequency. 
We  use  the  general  random-phase-approximation  (RPA) 
formula12  applied  to  our  2D  case  (<5p(x)  “*5p(r)5(r)l: 

Sp(r)  ~  f  dr'x^y)  t/e<l(r,)+  f  dr”-- .  (6) 

where  *(r,r')  is  the  susceptibility,  which,  for  circular  sym¬ 
metry,  has  the  form, 


*(r,r l  —  ’Egnj  £  - — — ^  , 

n.l  n‘.l+m  £n.l~en‘J  +  ni  +  fim  +  ty 


with  m  “l"  —  I  ,f„  t  the  electron  distribution  function,  and 
,  y  a  phenomenological  broadening  factor. 

Using  the  two-dimensional  Fourier  expansion  of 
I /I r'  —  r"|  and  expanding  Sp( r)  and  *(r,r')  as 

<5p(r)  —  YlSpni(r)ein‘° ,  (8) 


*(r,r  )“Z  y,„(r,r')e 


—  O') 


we  find 


5p, „(/•)-  Jo  dr'r'lxm(r,r')V^{(r‘)+8pm(r')pm(ry)] . 

(10) 

where 

F"'(r)  =  C*  deV'"(r)e""\  (11) 
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(13b) 


where  J,„(.v),  T(.t),;  F\(a,b,c:d)  are.  respectively,  the 
Bessel,  gamma,  and  Hypergeometric  functions.  The  hy- 
pergeometric  function  entering  Eq.  (13)  can  always  be  ex¬ 
pressed1’  as  a  linear  combination  of  the  complete  elliptic 
integrals  K(z )  ”  r  n-<F\ (  v ,  t  .  I;r 2)  and  E{z)mm'-n^ 
xF,(-l,  'r,Uz2i 

It  is  clear  that  the  solutions  of  the  integral  equation,  Eq. 
(10),  corresponding  to  different  m  are  independent. 
Furthermore,  only  those  modes  whose  w’s  correspond  to 
nonvanishing  Fourier  components  of  the  external  potential 
•  (Eq.  (1 1 )],  couple  to  the  external  radiation.  For  radiation 
whose  wavelength  is  much  lurger  than  the  dot  size,  the 
p  external  electric  field  can  be  considered  as  spatially  uni* 
r  form.  Choosing  this  field  to  be  polarized  along  the  x 


direction  so  that  Fe''(r)  —E^'x -E^VcosO.  Eq.  (11)  be¬ 
comes 

VS?(r) -E'"rn(SmJ  +Sm.  - . ) .  (14) 

Thus,  only  modes  with  w  ±  1  couple  to  the  external  ra¬ 
diation.  For  j5~0  these  modes  are  degenerate  (i.e.. 
8p\(r)  *5p_|(r)],  and  Eq.  (10)  need  only  be  solved  for 
wl.  For  this  case. 

-KlK(r‘/r"-E(r'/r")l  r'Kr"  (15a) 

R  i  (r*  r")  “ 

~\K{r"lr')-E{r"tr')V  r'>ru.  U5b> 
nr 

We  solve  Eq.  ( 1 0)  by  direct  matrix  inversion,  employing  a 


modified  Simpson’s  quadrature  method  to  account  for  the 
weak  singularities  of  the  elliptic  integral,  K  at  r'-r". 

From  the  induced  charge  density  we  calculate  the  cross 
section  for  photon  absorption,  aia),  given  by  a (a) 
■*4^Wr)Ima(ft)),  where  the  polarizability  a(a)  is 
defined  as  a(a)  ’m(\/Ecxl)fdrx8p(r).  In  the  long- 
wavelength  limit, 

drr2[8p\{r)  +  8p-\{r)] .  (16) 

As  an  example  we  consider  the  case  of  a  GaAs  dot 
structure  for  which  me  -0.0665m0  and  k  —  1 3. 1 .  We  take 
the  disk  radius,  /?  — 1000  A,  and  a  jellium  density, 
n+—4x|0n  cm-2,  which  corresponds  to  — - 1 25  positive 
charges  in  the  disk.  For  the  moment,  we  neglect  the  sur¬ 
face  states.  The  number  of  electrons,  ne,  is  varied  from  1 
to  30.  The  jellium  potential  V+(r ),  is  shown  in  Fig.  1 
(dotted  line).  Near  the  origin  this  potential  is  parabolic, 
as  can  be  verified  by  expanding  Eq.  (2)  for  small  r.  For 
larger  r  the  curvature  increases  above  that  of  a  parabola. 

The  solid  and  dashed  lines  in  Fig.  1  illustrate  the  self- 
consistent  ground-state  potential  and  corresponding  elec¬ 
tron  density  for  n,  ■“30.  Note  that  this  potential  is  flat  in 
the  region  of  large  electron  density  where  the  electrons 
screen  the  positive  background. 14  The  Fermi  energy, 
£/~6.3  meV,  15  spin-degenerate  levels  are  occupied  and 
the  energy  spacing  between  the  highest  occupied  and 
lowest  unoccupied  levels  is  around  1.5-2  meV. 

The  FIR  absorption  is  shown  in  Fig.  2  for  ne  “2  (dotted 
line),  n?“6  (dashed  line),  and  nf”30  (solid  line).  We 
take  the  broadening  factor,  y=0.1  meV.  For  all  cases, 
the  spectrum  is  dominated  by  a  single  peak.  For  small 
electron  number,  this  peak  is  very  close  to  the  frequency, 
ft»o,  of  the  parabolic  well  defined  by  the  curvature  near  the 
origin  of  the  jellium  potential,  and  it  is  essentially  in¬ 
dependent  of  n, ,  as  has  recently  been  observed  in  InSb  dot 
structures.2  Specifically,  expanding  Eq.  (2)  to  order  r 2 
yields,  o)6—e2Jtn+/KmeR.  In  this  case,  we  find  wo“4 
meV,  compared  to  <u— 4.2  meV  for  the  n^=*2  and  w,.  ”6 
absorption  peaks. 

We  note  that  the  absorption  frequency  is  significantly 
larger  than  the  energy-level  spacings  of  the  self-consistent 
ground-state  potential,  as  has  been  observed  in  expcri- 


0  200  400  600  800  1000  1200 
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FIG.  1.  Bare  quantum  dot  potential,  V  r,  (dotted  line),  as  de¬ 
scribed  In  the  text.  Also  shown  is  the  self-consistent  potential, 
F(r)  (solid  line)  and  electron  charge  density,  «(r)  (dashed  line) 
for  n,  "30. 
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FIG.  2.  FIR  response  of  the  dot,  cr(ft>),  in  the  absence  of 
magnetic  field  for  n,  —2, 6,  and  30. 


ment.2  This  large  depolarization  shift  cannot  be  assigned 
to  any  particular  single-particle  transition  since  many 
transitions  are  involved.  In  light  of  the  above  discussion, 
it  is  clearly  more  appropriate  to  associate  the  absorption 
frequency  with  the  energy  spacings  of  the  jellium  poten¬ 
tial. 

For  finite  magnetic  field,  the  orbital  degeneracy  is  lifted 
so  that  the  solutions  for  m“  +  l  and  m“  — I  are 
different.  This  leads  to  a  splitting  of  the  single  absorption 
peak  into  two  peaks,  which  we  identify  as  w+  and  «>-,  as 
shown  in  Fig.  3  for  the  case  of  ne  “6.  As  B  increases  the 
o)  —  mode  decreases  in  strength  while  the  tu+  mode  in¬ 
creases  in  strength.2  For  small  electron  number,  these 
peaks  occur  close  to  the  frequencies  appropriate  for  para¬ 
bolic  confinement: 

/ 

(0+  “  W+  j(l)})U2±  T  (Or  ,  (17) 

as  shown  in  the  inset  of  Fig.  3. 

For  large  electron  number,  the  absorption  shifts  to 
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FIG.  3.  FIR  response,  o(w),  as  a  function  of  magnetic  tieM. 
B.  Inset  shows  the  energies  of  the  absorption  peaks,  «*+  and  «*- 
(solid  circles)  vs  i»,  compared  to  the  strictly  parabolic  result 
(dashed  lines)  of  Eq.  (18)  with  woM  meV. 


higher  frequency  reaching  —15%  deviation  from  mo  for 
m,  "30  (solid  curve  of  Fig.  2).  A  similar  upward  shift  in 
absorption  frequencies  occurs  for  finite  magnetic  field. 
This  shift  is  related  to  the  changing  curvature  of  the  jelli- 
um  potential  with  r.  As  ne  increases  so  does  the  spatial 
extent  of  the  total  electron  density,  which  increasingly 
probes  the  nonparabolic  region  of  the  potential.  For  small 
deviations  from  parabolicity,  this  idea  can  be  given  quan¬ 
titative  formulation  by  taking  the  absorption  frequency, 
to  be  defined  by  the  average  curvature  of  the  jellium 
potential: 


2m,  R~  —  r~ 


with  the  average,  <•••),  taken  over  the  electron  charge 
density,  n(r).  Applying  this  formula  for  n,  “ 30,  we  find 
<ueff"4.6  meV,  in  excellent  agreement  with  the  RPA  re¬ 
sult15  (solid  line  in  Fig.  2). 

We  now  comment  on  the  effect  of  occupied  surface 
states.  It  is  evident  from  a  small  r  expansion  of  Eq.  (3) 
that  the  inclusion  of  surface  states  increases  the  curvature 
of  the  confining  potential  and,  as  a  result,  enhances  the 
resonance  frequency.  This  in  turn  enhances  slightly  the 
effects  of  nonparabolicity,  for  a  given  For  example,  for 
n,  “30  and  N  95,  we  find  <uo“4.7  meV  and  the  reso¬ 
nance  at  ro“5.5  meV,  an  upward  shift  from  (oq  of  17%. 

We  have  also  investigated  the  significance  of  ex¬ 
change-correlation  effects  by  including,  in  the  standard 
fashion,  an  exchange-correlation  model  potential, 
k\c ([«(/-)]),  in  the  local-density  approximation  (LDA) 
appropriate  for  a  two-dimensional  system.16  Specifically, 
for  the  calculation  of  the  ground  state,  Vxc  is  added  to  the 
yea  entering  Eq.  (4).  Similarly,  for  the  response  calcula¬ 
tion,  a  term  (dyK/dn)Sp( r')  enters  the  large  parentheses 
of  Eq.  (6).  We  find  a  small  downward  shift  of  the 
response  peak  whose  magnitude  decreases  with  increasing 


n,  and  is  at  most  10%  of  mo-  It  has  been  pointed  out,'6 
however,  that  two-dimensional  calculations  for  the 
exchange-correlation  energy  overestimate  substantially 
the  actual  exchange-correlation  effects  in  a  quasi-two- 
dimensional  system  with  finite  z  extent.  Therefore,  we 
conclude  that  these  effects  do  not  mask  the  nonparabolici¬ 
ty  effects  and  modify  only  slightly  the  RPA  results  ob¬ 
tained  above.  We  also  have  found  that  for  the  case  of 
strictly  parabolic  confinement  both  the  RPA  and  LDA 
calculations  produce  a  resonance  close  (within  1%  for 
n,  "30)  to  the  interlevel  separation  of  the  parabola. 

As  stated  above,  modes  with  differing  values  of  m  do 
not  couple  [see  Eq.  (10)]  because  of  the  circular  symme¬ 
try  of  the  dot  potential.  However,  once  this  symmetry  is 
broken,  modes  with  different  m  can  couple,  and  this  cou¬ 
pling  should  be  reflected  through  the  appearance  of  addi¬ 
tional  peaks  in  the  FIR  absorption  spectrum.  This  mecha¬ 
nism  provides  a  possible  explanation  for  the  mode  cou¬ 
pling  observed  in  recent  FIR  measurements  on  square- 
shaped  dots. 1 

In  summary,  we  have  developed  a  theory  to  describe  the 
electromagnetic  response  of  quantum  dots,  including  mag¬ 
netic  field.  We  find  that  for  small  electron  number  the 
FIR  absorption  spectrum  corresponds  to  that  associated 
with  parabolic  confinement.  For  large  electron  number, 
an  upward  shift  in  the  resonance  frequencies  occurs  as  the 
electron  density  probes  increasingly  the  larger  curvature 
of  the  dot  potential.  Furthermore,  additional  effects  of 
nonparabolic  confinement  such  as  multiple  mode  absorp¬ 
tion  and  mode  couplings  can  occur  if  the  dot  structures 
are  designed  with  strongly  reduced  symmetry  and  electron 
occupancy  that  is  sufficiently  high  to  probe  the  edge  asym¬ 
metries. 
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